Energy transfer among distant quantum systems in spatially shaped laser fields: 
Two H atoms with the internuclear separation of 5.29 nm (100 a.u.) 
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The quantum dynamics of two distant H atoms excited by ultrashort and spatially shaped 
laser pulses is studied by the numerical solution of the non-Born-Oppenheimer time-dependent 
Schrodinger equation within a three-dimensional (3D) model, including the internuclear distance R 
and the two z coordinates of the electrons, Z\ and z%. The two ID hydrogen atoms, A and B, are 
assumed to be initially in their ground states with a large (but otherwise arbitrary) internuclear 
separation of R = 100 a.u. (5.29 nm). Two types of a spatial envelope of a laser field linearly po- 
larized along the z-axis are considered: (i) a broad Gaussian envelope, such that atom A is excited 
by the laser field predominantly, and (ii) a narrow envelope, such that practically only atom A is 
excited by the laser field. With the laser carrier frequency u = 1.0 a.u. and the pulse duration 
t p = 5 fs, in both cases an efficient energy transfer from atom A to atom B has been found. The 
ionization of atom B achieved mostly after the end of the laser pulse is close to or even higher than 
that of atom A. It is shown that with a narrow spatial envelope of the laser field, the underlying 
mechanisms of the energy transfer from A to B and the ionization of B are the Coulomb attraction 
of the laser driven electron by the proton of atom B and a short-range Coulomb repulsion of the 
two electrons when their wave functions significantly overlap in the domain of atom B. In the case 
of a broad Gaussian spatial envelope of the laser field, the opposite process also occurs, but with 
smaller probability: the energy is transferred from the weakly excited atom B to atom A, and the 
ionization of atom A is also induced by the electron-electron repulsion in the domain of atom A due 
to a strong overlap of the electronic wave functions. 



I. INTRODUCTION 



The laser driven dynamics of distant quantum systems 
can depend in general on the spatial and temporal en- 
velope of the applied laser field. For example, at a gas 
pressure of 1 atm., the interparticle distance is about 
100 a.u. (5.29 nm). If such a system, e.g. composed 
of H atoms, is excited by a laser field with the carrier 
frequency uj = 1.0 a.u., corresponding to the ground- 
state energy of H2 at a large internuclear distance R, the 
wavelength is A = 861 a.u. (45.56 nm). If the laser field 
is focused within the diffraction limit onto a spot with 
the width A, the Gaussian spatial envelope of the field 
may result in quite different electric field strengths for 
H atoms separated by about 100 a.u., especially at the 
edges of the Gaussian spatial envelope. Although the 
H atoms are far away from each other, their electron- 
electron interaction should not be a priori neglected, es- 
pecially upon their excitation by the laser field, because 
the electronic wave functions extend and vanish, strictly 
speaking, only at infinity. Therefore, the energy transfer 
among distant quantum systems, similar to that studied 
in [1H8J, can be anticipated to occur in spatially shaped 
laser fields as well. For ultrashort laser pulses, containing 
only few optical cycles, one must also consider the carrier 
envelope phase (CEP) of the pulse [9[. 

The long-range energy transfer from an excited atom 
to its neighbor has been recently studied by Cederbaum 
et al. for molecular clusters [1] and is known as the 



interatomic Coulombic decay (ICD). Nowadays, ICD is 
well established also experimentally for inner-valence ex- 
citation of many electron systems I* 1 recent work 
Q, ICD was demonstrated experimentally for a helium 
dimer. Since helium atoms have no inner- valence elec- 
trons, a different type of ICD is operative for this case. 
It was thus concluded in [tJ that since ICD in a he- 
lium dimer takes place at interatomic distances up to 
Rj 12 a.u., no overlap of the electronic wave functions is 
required for the process. 

The present work is addressed to a quantum system 
composed of two H atoms with the initial internuclear 
separation of 100 a.u. (5.29 nm) which is excited by spa- 
tially shaped laser pulses: spatially broad pulses exciting 
both H atoms, and spatially narrow pulses exciting only 
one H atom of the entire H-H system. The relative sim- 
plicity of the H-H system under consideration (similar to 
that used in [l(|) makes it possible to treat the long-range 
electronic motion explicitly together with the nuclear mo- 
tion such as to reveal the role played by the electron- 
electron interaction and by the overlap of the electronic 
wave functions. An example of long-range laser-induced 
electron transfer (LIET) in the one electron linear H + - 
Hj atom-molecule system has been treated previously 
within the Born-Oppenheimer approximation [11] . Long- 
range charge and energy transfer can occur also in large 
molecular systems, as described recently in Ref. 12j and 
references therein. 

The following two types of H-H systems will be dis- 
tinguished in the present work: (i) a 'molecular' H-H 
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system, representing an elongated configuration of the 
H 2 molecule, similar to that studied recently in [l3[ for 
long-range entanglement, and (ii) an 'atomic' H-H sys- 
tem, representing two distant H atoms. Accordingly, the 
initial state of a molecular H-H system is assumed to 
be entangled by spin exchange and represented by the 
Heitler-London symmetric product of atomic wave func- 
tions, while the initial state of an atomic H-H system 
is not entangled - it is a direct-product state of atomic 
wave functions. In both cases the excitation of H-H is 
accomplished by laser pulses with (i) a broad Gaussian 
spatial envelope, such that both H atoms are excited by 
the laser field, with atom A being excited predominantly, 
and (ii) with a narrow spatial envelope, such that only 
atom A is excited by the laser field. 

The paper is organized as follows. The model of the 
H-H system and techniques used are described in Sec. II. 
Excitation, energy transfer, and ionization of an unen- 
tangled atomic H-H system are presented in Sec. III. Sec- 
tion IV is devoted to the laser-driven dynamics of an en- 
tangled molecular H-H system. The results obtained are 
summarized and discussed in the concluding Section V. 



II. MODEL, EQUATIONS OF MOTION, AND 
TECHNIQUES 

Within the 3D four-body model of H-H excited by 
the temporally and spatially shaped laser field the total 
Hamiltonian Ht is divided into two parts, 

H T (R,z 1 ,z 2 ,t) = H s (R,z 1 ,z 2 ) + H SF (z 1 ,z 2 ,t), (1) 

where Ha,(R,zi,z 2 ) represents the H-H system and 
Hsf(z\, z 2 ,t) describes the interaction of the system with 
the laser field. The applied laser field is assumed to be 
linearly polarized along the z-axis, the nuclear and the 
electronic motion are restricted to the polarization direc- 
tion of the laser electric field. Accordingly, two z coor- 
dinates of electrons, z\ and z 2 , measured with respect to 
the nuclear center of mass, are treated explicitly together 
with the internuclear distance R. A similar model has 
been used previously in 10] for the H 2 molecule, where 
each particle, electron or proton, is treated in ID, i.e., z 
and R. 

The total non-Born-Oppcnhcimer system Hamiltonian 
(employing a. u.: e = fi = m e = 1) reads 



H S (R, zi,z 2 ) 



1 d 2 
m p dR 2 



V PP (R) 



polarization terms are neglected. The Coulomb poten- 
tials in Eq. @ read 

V PP (R) = V ee {zi,z 2 ) = , 1 =5= =, 

R ^{zx- z 2 ) 2 +a 



V ep (z k ,R) = -- 



y/(z k - R/2) 2 + f3 y/(z k + R/2f+p 

(3) 

where k = 1,2, and the regularization parameters, a — 
0.1 x 10 -3 and /3 = 1.995, have been chosen (similar 
to previous work [13]) such as to reproduce the ground- 
state (GS) energy of the H-H system at R = 100 a.u. 
(EgS H = _1. a .u.). 

The interaction of the H-H system with the laser field 
is treated within the semiclassical electric dipole approx- 
imation by the Hamiltonian 



#sf (zi , z 2 , t) = - - (1 + 7 ) 22 F ( z k)zk , 



(4) 



fe=i 



where 7 = (1 + 2m p ) , A(t) is the vector potential, c is 
the speed of light, and F(z) is the spatial-shape function, 
or envelope, of the laser field. In the general case we 
set F(zi) — F(z 2 ), except for some model simulations 
detailed in Sec. IV. 

The vector potential, A(i), is chosen in the following 
form: 



A(t) = — £q sin 2 (nt/tp) cos(ujt + . 



(5) 



where £q is the amplitude, t p is the pulse duration at the 
base, to is the laser carrier frequency, and <f> is the carrier- 
envelopephase (CEP). Note that it has been shown pre- 
viously Q that the carrier phase of the laser pulse is 
important only for pulses having less than 15 optical cy- 
cles. The definition of the system-field interaction by 
Eq. (UJ) via the vector potential, suggested in [l3j], as- 
sures that the electric field £(t) = —~dA(t)/dt has a 

vanishing direct-current component, J* p £(t)dt — 0, and 
satisfies Maxwell's equations in the propagation region. 

It is suitable to define, on the basis of Eqs. Q and ([5]), 
the local effective-field amplitudes for atoms A and B as 
follows: 



£ A = £ F(z A ), £ B = £ F(z B ), 



(6) 



where za = —50 a.u., zb = 50 a.u.. The respective time- 
dependent electric fields acting on atoms A and B read 



.-A,B/ 



(t) = £g' B [sin 2 {irt/tp) sin(wi + ■ 



fe=i 



2/i c dz 2 



2 + Vep(z k ,R) 



where m p is the proton mass, fi c 
is the reduced electron mass, and non-diagonal mass- 



+ V ee (zi,z 2 ), 



2m p / (2m p 



(2) 



1) 



cut. 



■ sin(27rf/i p ) cos(wi + <p)]. 



(7) 



The first term in Eq. ([7]) corresponds to a laser pulse with 
a sin 2 -type temporal envelope, while the second, the so- 
called 'switc hing ' term, appears due to the finite pulse 
duration E3, Km. 
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The 3D time-dependent Schrodinger equation for ID 
electrons (z±, z 2 ) and ID protons (R), 



d 

i — ty = [H s {R,z u z 2 ) 



(8) 



has been solved numerically with the propagation tech- 
nique adapted from [loT ] for both electron and proton 
quantum motion. In particular, calculations for the elec- 
tron motion have been performed by using 200-point non- 
equidistant grids for the Hermite polynomials and corre- 
sponding weights for the numerical integration on the in- 
terval (—00, 00) for the Z\ and z 2 coordinates. For the nu- 
clear coordinate R, a 256-point equidistant grid has been 
used on the interval [75 a.u., 125 a.u.]. The time-step of 
the propagation was At = 0.021 a.u. (1 a.u. —24 asec). 

The wave functions of the initial states have been 
obtained by numerical propagation of the equation of 
motion (|8]) in imaginary time without the laser field 
(£o=0). 

Upon excitation of the H-H system by the laser field, 
the electronic wave functions of its atomic A and B parts 
may overlap. In order to study the energy transfer from 
A to B, the respective 'atomic' energies, E^(i) and Eft(t), 
are defined on the basis of Eqs. © and © such that their 
sum always gives the correct total energy of the entire 
H-H system. The definitions of 'atomic' energies E^(t) 
and Eb (t) for the H-H systems of atomic and molecular 
origin are different due to the different symmetry and 
entanglement of the respective wave functions and will 
be specified in the following sections. 

The ionization probabilities for atoms A and B have 
been calculated from the time- and space-integrated out- 
going fluxes separately for the positive and the negative 
directions of the Z\ and z 2 axes at Z\ 2 = ±91 a.u.. 
Specifically, we calculated four ionization probabilities: 
Ia(^i = —91 a.u.) and Ia(^2 = —91 a.u.) for atom A, 
and Ib(zi = 91 a.u.) and Ib(z 2 = 91 a.u.) for atom B. 
At the outer limits of the z-grids, absorbing boundaries 
have been provided by imaginary smooth optical poten- 
tials adapted from that designed in Similar optical 
potentials have been also provided for the i?-axis but, 
in practice, the wave-packet never approached the outer 
limits of the i?-grid. 



III. EXCITATION OF H-H FROM AN 
UNENTANGLED DIRECT-PRODUCT INITIAL 
STATE 

The spatial part of the initial unentangled direct- 
product ground-state wave function of H-H of an atomic 
configuration used in the imaginary time propagation is 
defined with the unsymmetrized Heitler-London electron 
wave function as follows: 



9(R, z u z 2 ,t = 0) 

where ®is A {zi) = 
*is B (^) = e-^- z ^ 
a proton Gaussian function centered at R = 100 a.u.. 



= *lS A (*0*lS B (*2)*G(fl), (9) 

e -\ z i- z *\ at z A = -50 a.u., 
at zb = 50 a.u., and 'I'g(-R) is 



The imaginary-time propagations have been performed 
with a reduced, 'non-interacting' atoms, version of the 
system Hamiltonian @ wherein the nuclear V PP (R) and 
electronic V ee (zi, z 2 ) terms were omitted and only the 
Coulombic interaction V e p(zu2), R) of each electron with 
its nearest proton was taken into account along with the 
three kinetic-energy terms of Eq. ©. The results ob- 
tained gave a direct-product initial state. For the sake of 
comparison, we have also performed the imaginary-time 
propagations with the complete system Hamiltonian ^ . 
The results obtained proved to be practically identical 
to those obtained in the case of two non-interacting H 
atoms. The energy difference, for example, is less than 
1.8 x 10~ 5 a.u., implying that the Coulombic interaction 
V ee (z\, z 2 ) of two H atoms in their ground states is neg- 
ligible at R = 100 a.u.. 
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FIG. 1. The initial unentangled direct-product state of the 
H-H system and the spatial envelopes of the applied laser 
pulses, (a) - electron probabilities: -P(zi) is the probability 
to find electron ei initially belonging to atom A at z = Z\ with 
any values of the other two coordinates, R and Z2\ probability 
P{z2) has a similar meaning for electron e2 initially belonging 
to atom B; (b) - spatial envelopes of applied laser fields: the 
broad Gaussian envelope of Eq. (|14l) and the narrow envelope 
of Eq. Ip)). 



The initial unentangled direct-product state of the H-H 
system representing two non-interacting H atoms, A and 
B, is presented in Fig. 1(a) by the electron probabilities 
P(zi) and P(z 2 ), which are defined as follows: 



P(z 1 )= J dR J dz 2 \y(R,z l7 z 2 )\ 2 



(10) 



for electron e\ initially belonging to atom A with proton 
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Pa, and similarly, 

P(z 2 ) = J dpj dz 1 \^{R,z 1 ,z 2 )\ 2 , (11) 

for electron e 2 initially belonging to atom B with proton 
Pb ■ These electron probabilities give the overall probabil- 
ity to find an electron at a specified point of the z-axis at 
any position of the other electron and at any internuclear 
distance. 

After calculation of the unentangled initial state of the 
atomic H-H system, its laser-driven quantum dynamics 
in real time has been explored with the complete system 
Hamiltonian of Eq. ([2]). The spatial envelopes of the 
applied laser pulses, Gaussian and narrow, defined by 
Eqs. ([TJ]) and Eq. (fTB")) below, are presented in Fig. 1(b) 
to illustrate the local effective-field strengths acting on 
atoms A and B. 

Two possible choices for the laser carrier frequency are: 
uj = 0.5 a.u. (corresponding to the ground-state energy 
of an H atom), and uj = 1.0 a.u. (corresponding to the 
ground-state energy of H-H at large R). Our numerical 
simulations have shown that at uj — 0.5 a.u., an efficient 
ionization of A takes place, while the energy transfer from 
A to B and the ionization of B are not efficient. In con- 
trast, at ui = 1.0 a.u., the energy transfer from A to B 
and the ionization of B are efficient. This implies the ex- 
istence of an optimal laser frequency suitable for the most 
efficient energy transfer. This issue will be addressed in 
a forthcoming work; below we present the results ob- 
tained at uj = 1.0 a.u.. With the laser pulse duration 
t p = 5 fs, as used throughout this work, the number of 
optical cycles during the pulse, N c — ujtp/2n, is about 
33. Therefore, the carrier phase of the laser field </> is not 
important [9| and set equal to zero in our simulations 
below {(j) = in Eqs. © and ©]. 

In order to study the energy transfer from A to B, the 
respective 'atomic' energies, E A (t) and Es(t), have been 
defined on the basis of Eqs. ([2]) and (J3]) as follows: 



E A (t) = (*(*) 



_1 

2m p dR 2 



+ \v pv {R) + \v cc {z u z 2 ) 



2^ c dz\ 



E 



\ y/{z k + R/2)*+0 



*(*) 



(12) 



for atom A and similarly, 

1 d 2 



E B (t) = (*(*) 



2m, dR- ' 2 V ^ + 2 V "^ 



i d 2 A i 

2Me dz 2 ^ ^/( Zk - R/2Y + (3 



*(*)>, (13) 



of the proton-proton and the electron-electron interac- 
tion are assumed to be equally shared between atoms A 
and B. The kinetic energy of electron e\ and the energy of 
Coulombic interaction of both electrons, e\ and e 2 , with 
proton pa are entirely assigned to atom A. Similarly, the 
kinetic energy of electron e 2 and the energy of Coulombic 
interaction of both electrons with proton pb are entirely 
assigned to atom B. The sum of 'atomic' energies always 
gives the correct total energy of the entire H-H system. 
Notice, in particular, that the fact that the kinetic energy 
of electron ei/e 2 is entirely assigned to atom A/B corre- 
sponds to the initial electron probabilities P(z\)/ P(z 2 ) 
in the unentangled atomic state [Fig. 1(a)]. Indeed, elec- 
tron ei is localized in the vicinity of proton p A of atom A 
and electron e 2 is localized in the vicinity of proton p-Q of 
atom B. Therefore, if e.g. the laser pulse with a narrow 
spatial shape [Fig. 1(b)] is used to excite the extended H- 
H system, only the 'atomic' energy E A (t) will increase, 
while E-ait") will not be affected. 

The left panel of Fig. 2 presents the dynamics of the 
atomic H-H system excited by the laser pulse with the 
Gaussian spatial envelope 



F G (z) = cxp{-[(z-z )/A] 2 }, 



(14) 



where A = 861 a.u. (45.56 nm) and the laser field is 
assumed to be focused onto a spot centered at zq, where 
z Q = -1.5A [z =-1291.5 a.u. (-68.34 nm)]: cf. Fig. 1(b). 
The laser pulse parameters are: uj = 1.0 a.u., t p = 5 fs, 
and So — 1.0 a.u. (£ = 5.14 x 10 9 V/cm, intensity 1$ = 
3.5 x 10 16 W/cm 2 ). Accordingly, the spatial envelope 
function F(z) in the interaction Hamiltonian of Eq. (|4| 
is set equal to Fq(z) of Eq. (fl4|) . and the effective-field 
amplitudes for atoms A and B defined by Eq. ^ are: 
£fr = 0.125 a.u. (intensity l£ = 5.48 x 10 14 W/cm 2 ) and 
£jf = 0.088 a.u., (i? = 2.72 x 10 14 W/cm 2 ), implying a 
dominant excitation of atom A, which is clearly observed 
in Fig. 2(a). 

From Fig. 2(a) we see that the 'atomic' energy Exit) 
is controlled by the applied laser pulse: it increases in 
the first half of the pulse and decreases to the end of the 
pulse. In contrast, the 'atomic' energy E-B,{t) does not 
follow the applied laser pulse: it slowly increases up to 
the end of the laser pulse and even exceeds E^t) shortly 
after the end of the pulse. A similar behavior is found for 
the 'atomic' ionizations Ia and Ib presented in Fig. 2(b): 
while the laser-induced ionization Ia rises fast in the sec- 
ond half of the pulse, the ionization probability Ib sharply 
rises after the end of the pulse. Such a behavior will be 
referred to below as a 'sequential' ionization, implying a 
delayed interaction due to LIET. The time-dependent ex- 
pectation values (R(t)), {z\(t)) and (z 2 (t)) are presented 
in Fig. 2(c). The time-dependent deviations of the ex- 
pectation values (zi(t)} and (z 2 (t)), 

A«i(t) = (z 1 (t))-{z 1 (t = 0)), 



for atom B. From Eqs. (|T2|) and ([13]) we see that the ki- 
netic energy of nuclear motion and the potential energies 



Az 2 (i) = (32(f)) - (z 2 (t = 0)), 



(15) 
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FIG. 2. Quantum dynamics of H-H excited from an unentangled direct-product initial state by spatially shaped laser fields: a 
Gaussian spatial envelope (a-d) and a narrow spatial envelope (e-f). (a) and (e) - the energy transfer; (b) and (f) - 'atomic' 
ionizations [ionization probabilities Ia(zi = 91 a.u.) and Ib(z2 = —91 a.u.) are of the order of 10 -6 and are not shown] ; (c) and 
(g) - expectation values {R), (z\) and (22) (deviations of the expectation values, Azi and Az2, defined by Eqs. (|15|) are scaled 
up by a factor of two); (d) and (h) - local effective laser fields S A (t) [Eqs. ([6| and l[7|)] acting on atom A. 



plotted in Fig. 2(c) are scaled up by a factor of two. We 
see from Fig. 2(c) that Azi(i) > 0, implying attraction of 
electron e\ by proton pe- In contrast, Az%(t) < 0, imply- 
ing attraction of electron e2 by proton It is also seen 
from Figs. 2(a), 2(b) and 2(c) that both the energy trans- 
fer from A to B and the 'sequential' ionization of B in the 
positive direction of the Z2-axis correlate with a small de- 
crease of the spatial separation of electrons, \(z2) — { z i)\, 
such that the electron-electron repulsion (EER) becomes 
effective. 



In order to clarify these findings further, we have per- 
formed a similar numerical simulation with a narrow spa- 
tial envelope of the laser pulse acting practically only on 
atom A of the entire H-H system. Specifically, the fol- 
lowing sin 2 -type spatial envelope of the laser pulse was 
used for z = Z\ and z = z 2 : 



Fs{z) 



sm 



7r(z - Z a ) 



Zb 



Za < Z < Z b , 



(16) 



where z„ 



-60 a.u., z — —40 a.u., and F^(z) = oth- 
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erwise, as illustrated in Fig. 1(b). Accordingly, the spa- 
tial envelope function F(z) in Eq. (Q| is set equal to F N (z) 
of Eq. (fT6|) . The amplitude of the pulse, £q = 0.02 a.u. 
(Jo = 1-4 x 10 13 W/cm 2 ), is chosen such that the 'atomic' 
energy E\(t) gained in the field with the narrow spa- 
tial envelope [Fig. 2(e)] is close to that gained with the 
Gaussian spatial envelope [Fig. 2(a)]. Note that, defined 
by Eq. (O, the effective-field strength £ A (t) of the laser 
pulse with the narrow spatial envelope [Figs. 1(b) and 
2(h)] is much smaller than that of the laser pulse with 
the Gaussian spatial envelope [Figs. 1(b) and 2(d)], but 
the energy transfer from A to B is by about 30% more 
efficient [Fig. 2(e)]. Also note that, due to the effective- 
field amplitude acting on atom B in the current case of 
the narrow spatial envelope of the pulse is £ g 3 =0, the ef- 
ficient energy transfer to atom B is attributed entirely to 
the electron-electron repulsion T4 e (zi,z 2 ) and electron- 
proton attraction 14 p (0 1 ( 2 ), R) in the vicinity of atom B 
( z i(2) = z b)' This also applies to the ionization of atom B 
in the positive direction of the z 2 -axis [Fig. 2(f)], which is 
not only by about 30% more efficient than that induced 
by the laser pulse with the broad Gaussian spatial enve- 
lope [Fig. 2(b)], but Ib is even larger than T\ at t > 7 fs 
[Fig. 2(f)], implying that the 'sequential' ionization Ib, 
induced by EER, is more efficient than the laser-induced 
ionization Ia- 

Taking into account that in the unentangled direct- 
product state electrons ei and e 2 are well localized 
on the z-axis [Fig. 1(a)], one can assume that the time- 
dependent expectation values (z\(t)) and (z 2 (i)} repre- 
sent the positions of electrons. A closer look at the data 
presented in Figs. 2(c) and 2(g) for {z\(t)) and (z2(t)) 
shows that the two electrons of H-H do not approach 
each other closer than 98 a.u. (5.18 nm). Therefore, at 
a first glance, the energy transfer from A to B and the 
'sequential' ionization of B take place, similarly to ICD 
[lH8j], without any noticeable overlap of the respective 
electronic wave functions. Nevertheless, a closer look at 
spatial distributions of the electronic wave functions illus- 
trated in Fig. 3 reveals that EER is in fact a short-range 
process which takes place in the vicinity of atom B at 
zb = 50 a.u.. Here, the Coulombic attraction of electron 
ei by proton pb effectively localizes electron e±. Accord- 
ingly, the overlap of the electronic wave functions has a 
sharp narrow maximum at zb = 50 a.u. (Fig. 3), which 
results in a strong EER. 

Electron probabilities P(zi) and P(z 2 ) presented in 
Fig. 3 correspond to the end of the 5 fs laser pulse with 
the narrow spatial envelope defined by Eq. (fT6|) . This 
pulse excites practically only electron e\ in the vicinity 
of za = —50 a.u., therefore the probability distribution 
P{z\) at t — 5 fs is broadened in comparison to the ini- 
tial one [Fig. 1(a)] in both Z\ < za and Z\ > za di- 
rections. The laser- induced extension of P(z\) into the 
domain of z\ < za gives rise to the ionization Ia on the 
negative direction of the z^axis. At z\ < —91 a.u., the 
wave function is absorbed by the imaginary optical po- 
tential. In contrast, at z\ > za, the laser-driven electron 
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FIG. 3. Electron probabilities P{zi) (solid line) and P(z2) 
(dashed line) at the end of the 5 fs laser pulse with a narrow 
spatial envelope [Eq. (|16[) and Fig. 1(b)] which excites only 
electron ei (coordinate Zi) in the domain of atom A. 



e\ reaches the domain z\ > where it is attracted and 
accelerated by proton pb, localized at zb = 50 a.u.. Due 
to the electron-proton Coulombic attraction, the prob- 
ability P(z\) has its local maximum at Z\ = zb, cor- 
responding to LIET of Ref. [TTJ] . Simultaneously, when 
electron e\ approaches the domain of zb = 50 a.u., where 
the initial probability -P(z 2 ) of electron e 2 has the global 
maximum [Fig. 1(a)], the electron-electron Coulombic re- 
pulsion becomes very strong. Therefore, the probability 
P(z 2 ) is extended into the domain of z 2 > zb (Fig. 3), 
giving rise to the ionization Ib in the positive direction of 
the z 2 -axis. Finally, the wave function is absorbed by the 
imaginary optical potential at z 2 > 91 a.u.. Note that 
in the current case of a narrow spatial envelope of the 
laser field [Fig. 1(b)] one can clearly distinguish between 
the laser-induced ionization on the negative direction of 
the zi-axis, represented in Fig. 2(f) by \a{z\ = —91 a.u.), 
and the 'sequential' EER-induced ionization on the pos- 
itive direction of the z 2 -axis, represented in Fig. 2(f) by 
Ib(^2 = 91 a.u.). Also note that the ionization probabili- 
ties \a{z\ = 91 a.u.) and Ib(^2 = —91 a.u.) are less than 
10" 6 at t = 5 fs. 
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FIG. 4. Electron probabilities P(zi) (solid line) and P(zi) 
(dashed line) at 2.5 fs after the end of the 5 fs laser pulse with 
a Gaussian spatial envelope which excites both electrons, with 
the excitation of electron ei (coordinate z%) being dominant 
in comparison to that of electron (coordinate 22). 
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In the case of the broad Gaussian spatial envelope of 
the laser pulse centered at z = — 1291. 5 a. u. (-68.34 nm), 
both electrons, e\ and e%, are excited by the laser field 
and therefore ionization in both positive and negative di- 
rections of the z-axis is induced by both the laser field and 
subsequently by EER. The electron probabilities P{z\) 
(solid line) and P{z2) (dashed line) for the Gaussian spa- 
tial envelope of the 5 fs laser pulse are plotted in Fig. 4 
at t = 7.5 fs, when both ionization probabilities plotted 
in Fig. 2(b) approach their maximum values. Although 
both electrons are excited by the laser field, the excitation 
of electron ei in the domain of atom A is about 1.4 times 
stronger than the excitation of electron ei in the domain 
of atom B. Accordingly, the local maximum of P(zi) at 
zb = 50 a.u. (Fig. 4) is almost 4 times higher than the 
local maximum of P{z2) at za = —50 a.u. and therefore 
the EER-induced ionization in the positive direction of 
the Z2-axis is stronger than the EER-induced ionization 
in the negative direction of the zi-axis. In contrast, the 
laser-induced ionization in the negative direction of the 
zi-axis is stronger than that in the positive direction of 
the Z2-axis. Indeed, it is seen from Fig. 2(b) that at the 
end of the 5 fs laser pulse, the ionization probability Ib 
is more than twice smaller than Ia, while Ib ~ Ia at 
t >10 fs. 

It is seen from Fig. 4 that both local maxima, P(zi) 
at z\ = —50 a.u. and P{z2) at Z2 — 50 a.u., are sig- 
nificantly smaller than the local maximum of P(z\) at 
z\ = —50 a.u. produced by the laser pulse with a narrow 
spatial envelope (see Fig. 3). Therefore the EER-induced 
ionization in the case of the broad Gaussian spatial en- 
velope of the laser pulse [Fig. 2(b)] is smaller than that 
in the case of the narrow spatial envelope of the pulse 
[Fig. 2(f)]. 

It can be concluded from Fig. 4 that in the general 
case, when both distant atoms are excited by the laser 
field, energy is transferred from A to B and from B to 
A, and EER-induced ionization occurs in both A and B 
parts of H-H along with the laser-induced ionization. If, 
for example, the laser pulse has a wide Gaussian spatial 
envelope and/or is centered at z = 0, the mutual energy 
transfers and the EER-induced ionization probabilities 
are substantial for both distant atoms A and B even at 
a large internuclear separation. 

Finally we note that the appearance of sharp local 
maxima of P{z\) at z\ = 50 a.u. and of P{z2) at 
Z2 = —50 a.u. (see Figs. 3 and 4) confirm spreading, 
derealization, and non-factorization of the initially fac- 
torized and well localized on the z-axis wave function 



of H-H [Fig. 1(a)] and can therefore [171 be treated as 
the emergence of long-range entanglement, or quantum 
non-local connection, at R = 100 a.u.. Recently this 
long-range entanglement attracted considerable interest 
both in theory and experiment [HI, [l8l - l2~i| . 



IV. EXCITATION OF H-H FROM AN 
ENTANGLED INITIAL STATE 

The H-H molecular system represents an elongated 
configuration of the H2 molecule. Therefore, its ini- 
tial electronic ground state is entangled via exchange 
@ E3, Hil- The spatial part of the initial entangled 
ground-state wave function of H-H for a singlet electronic 
state is given by 

%(R, Zl ,z 2 ,t = Q) = [*1S A (Z1)*1S B (2 2 ) 



(17) 



where $is A B (21,2) are defined by Eq. © and SS?g(R) is 
a proton Gaussian function centered at R = 100 a.u.. 
Imaginary time propagations have been performed with 
the complete system Hamiltonian ([2]). In the entangled 
initial state of the H-H system, electron probabilities 
P(zi) and P{z2) are identical to each other. 

For the simulations with this entangled initial state 
(fl"7|) . the 'atomic' energies E^(t) and Esit) are defined 
on the basis of Eqs. @ and ([3]) as follows: 



E A (t) = ( *(t) 



2m p dR 2 



fe=i 



1 1 



*(*)> (18) 



for atom A, and 
E B {t) = /*(*) 



1 d 2 



2m p dR 2 2 ^ 2fi c dz\ 2 



1 2 

-V ee (z u z 2 ) - 



£j ^{z k -R/2) 2 +p 



*(*)> (19) 



for atom B. From Eqs. (|18|) and (|19|) we see that the 
kinetic energies of nuclear and electronic motion as well 
as the potential energies of the proton-proton and the 
electron-electron Coulombic interaction are assumed to 
be equally shared between atoms A and B. The energy 
of Coulombic interaction of both electrons with proton 
Pa is assigned to atom A, while the energy of Coulombic 
interaction of both electrons with proton p& is assigned 
to atom B. The sum of these 'atomic' energies always 
gives the correct total energy of the entire H-H system. 
The choice of the electron kinetic energies in Eqs. (fT8|) 
and (TiTH) corresponds to the initial electron probabilities 
P{z\) and P(z2) in the entangled molecular state. Both 
electrons e\ and e 2 are localized with the 50% probability 
in the vicinity of proton p A of atom A and if the extended 
H-H system is excited e.g. by the narrowly shaped laser 
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FIG. 5. Quantum dynamics of H-H excited from the entangled initial state (|17[) by spatially shaped laser fields [Fig. 1(b)]: a 
broad Gaussian spatial envelope (left panel) and a narrow spatial envelope (right panel), (a) and (d) - the energy transfer; (b) 



and (e) - 'atomic' ionizations; [the total ionization probabilities I ^ 
effective laser fields £ A (i) [Eqs. ((6]) and ©J acting on atom A. 



and l£ otal 



are defined by Eq. ([20]) ]: (c) and (f) - local 



pulse of Fig. 1(b), both electrons give rise to the 'atomic' 
energy £?a(£), eacn with 50% probability. 

The quantum dynamics of H-H excited from the en- 
tangled initial state (|T7|) by spatially shaped laser pulses 
is presented in Fig. 5. The left panel corresponds to 
the broad Gaussian spatial envelope centered at zq — 
— 1291.5 a.u. (-68.34 nm) and the right panel corresponds 
to the narrow spatial envelope centered at za = —50 a.u. 
(-2.65 nm), see Fig. 1(b). For the sake of comparison, the 
laser fields used are the same as in the previous case of the 
unentangled direct-product initial state [see Figs. 1(b) 
and 2] . Specifically, the Gaussian spatial envelope is de- 
fined by Eq. (fT4")l , the narrow spatial envelope is defined 
by Eq. (|16p , the carrier frequency of the laser pulse with 
the sin 2 -type temporal envelope of Eq. (JJJ) is u = 1.0 a.u., 
and the pulse duration at the base is t p = 5 fs. 

From the comparison of the quantum dynamics of the 
unentangled state [Figs. 2(a) and 2(e)] to that of the 
entangled state [Figs. 5(a) and 5(d)] the following obser- 
vations are made. 

(i) The overall energy AEb transferred on a long 



timescale of t = 20 fs from atom A to atom B in 
the entangled state is very similar to that transferred 
from A to B in the unentangled direct-product state: 
A£b ~ 0.01 a.u. for the broad Gaussian spatial en- 
velope, and A£b ~ 0.013 a.u. for the narrow spatial 
envelope. 

(ii) The maximum 'atomic' energy Exit) gained during 
the laser pulse by atom A in the entangled state is smaller 
than that gained in the unentangled direct-product state. 
In contrast, the maximum 'atomic' energy E^{t) gained 
during the pulse by atom B in the entangled state is 
substantially larger than that gained in the unentangled 
direct-product state. 

(hi) Moreover, in the case of the entangled initial state, 
the 'atomic' energy Es,{t) is controlled by the laser pulse 
similarly to E^{t), even when only atom A is excited 
by the laser pulse with a narrow spatial envelope: the 
energy Eb (t) increases in the first half of the laser pulse 
and decreases at the end of the pulse, in contrast to the 
case of the unentangled direct-product initial state. This 
observation implies that the changes made by the applied 
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laser field to the entangled wave function in the domain 
of atom A at z 12 ~ —50 a.u. result in simultaneous 
changes in atom B at z\ t % ~ 50 a.u. due to the symmetry 
of the wave function by exchange. Such an entangled 
behaviour is very important for a long-range quantum 
communication among distant quantum systems jl9j |. 

(iv) After the end of the laser pulse, the 'atomic' en- 
ergies Ep^{t) and E&(t) demonstrate out-of-phase oscil- 
lations: slow oscillations in the case of the unentangled 
direct-product initial state [Figs. 2(a) and 2(e)] and slow 
oscillations modulated with very fast ones in the case of 
the entangled initial state [Figs. 5(a) and 5(d)]. The am- 
plitudes of both slow and fast oscillations of Exit) and 
E-g(t) decrease with time, which may indicate the for- 
mation of a quasi-stable configuration of the excited H-H 
system. 

It is instructive, before proceeding, to consider the 
laser-driven dynamics of the entangled H-H system in 
the initial stage of its excitation by the laser pulse with 
a narrow spatial envelope (|16[) which excites electrons 
ex and e 2 only in the domain of atom A. In Fig. 6, the 
laser-driven dynamics of entangled H-H is presented on 
the timescale of 1 fs. Time-dependent 'atomic' energies, 
Ex(t) and Esit), and expectation values of electronic 
coordinates, (z\(t)) and {z 2 {t)), are shown in Figs. 6(a) 
and 6(b) respectively. The time-dependent laser field is 
shown in Fig. 6(c). 

It is seen from Fig. 6(a) that the 'atomic' energies 
E\(t) and Es(t) oscillate out-of-phase with respect to 
each other, with Ex(t) being in-phase and E B (t) being 
out-of phase with the applied laser field [Fig. 6(c)]. 

As seen from Figs. 6(b) and 6(c) the electrons also 
follow the applied laser field. Note that expectation val- 
ues (zi(t)) and (z2(t)} in the entangled H-H system are 
identical and not distinguishable in Fig. 6(b). A perfect 
electron-field following at the laser carrier frequency be- 
ing as high as 1 a.u. is very interesting. Previously, the 
electron-field following on the level of expectation values 
of electronic coordinates have been explored only in the 
infrared [l5| and near-infrared [22| domains of the laser 
carrier frequency. The electron-field following at high 
laser frequencies is reminiscent of the well-known recolli- 
sion model of Corkum [23| . However, an important fea- 
ture of Fig. 6(b) is that electrons follow the applied laser 
field in-phase, while according to the theoretical model 
used in |23j and numerical results of [l5[ electrons follow 
the field out-of-phase: (z(t)) decreases when electric-field 
strength £{t) increases. A detailed study of the electron 
dynamics showed that, probably due to their finite al- 
beit very small mass, the electron do not react to the 
first half-cycle of the applied field at lu = 1 a.u. and 
follows the field in-phase at t > 0.15 fs. On the other 
hand, the out-of-phase electron-field following take place 
at lu < 0.1 a.u.. Similar results have been obtained for 
the unentangled H-H system excited by the laser pulse 
with a narrow spatial envelope: the electron e\ of atom A 
follows the applied field in-phase at lu = 1 a.u. and out- 
of-phase at lu < 0.1 a.u.. These results can be explained 
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FIG. 6. The initial (1 fs) stage of excitation of H-H from 
the entangled initial state [Eq. (I17|) ] by the laser pulse with 
a narrow spatial envelope [Eq. (|16p ] acting on atom A. (a) - 
time-dependent 'atomic' energies E\(t) (solid line) and E&(t) 
(dashed line); (b) - electron- field in-phase following; (c) - the 
laser field acting on atom A. 



as follows. At lu < 0.1 a.u., the electron is in the ground 
state, the polarizability is negative and, therefore, the 
electron follows the applied laser field out-of-phase. In 
contrast, at ui = 1 a.u., the electron is well above the ex- 
cited state, the polarizability changes the sign, and elec- 
tron follows the field in-phase. Similar behaviour has 
been observed in a recent work [2~i| for molecular ion H^J~ 
excited at wavelength A =800 nm (lu = 0.057 a.u.). At 
fixed R = 2 a.u., the electron is in the ground state and 
follows the field out-of-phase, whereas at R = 7 a.u., the 
laser carrier frequency lu — 0.057 a.u. is larger than the 
energy difference E\ sau — Ei sag , and the electron follows 
the field in-phase. 

Coming back to the laser-driven dynamics of the entan- 
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gled H-H system on the long timescale of 20 fs (Fig. 5), we 
note that ionization of atom B starts only in the second 
half of the laser pulse [see Figs. 5(b) and 5(e)]. Due to the 
symmetry of the entangled wave function, the ionization 
probabilities of electrons e\ and e\ are identical to each 
other both for the positive and the negative directions 
of the z-axes: Ia(^i = — 91 a.u.) = \a{z 2 = — 91 a.u.) 
and Ib(zi = 91 a.u.) = Ib(z2 = 91 a.u.). Therefore, in 
Figs. 5(b) and 5(e) the total ionization probabilities, 

T Totai _ Ia ( Zi = _ gi a u ) + Ia ( Z2 = _ gl a u ^ 

Igotal = j b ( Zi = gl a u -j + j b ( Z2 = 91 a u ^ ( 20 ) 

are plotted. From the comparison of Figs. 2(b) and 2(f) 
to Figs. 5(b) and 5(e) one can easily see that the total 
time-dependent ionization probabilities I^ otal and I^otai 
used in the case of the entangled initial state are almost 
equal (but not identical) to the time-dependent ioniza- 
tion probabilities Ia{z\ = —91 a.u.) and Ib(z2 = 91 a.u.) 
used in the case of the unentangled direct-product initial 
state. One can conclude therefore that the entanglement 
of the initial state of H-H, including its exchange symme- 
try, does not change its ionization probability in compari- 
son to the non-symmetric unentangled direct-product ini- 
tial state. 

The other consequence of the exchange symmetry of 
the entangled wave function is a very small spatial sep- 
aration of electrons all over the z-grid. The time- 
dependent expectation values (zi(t)) and (z2(t)} are al- 
most identical, both being close to on the timescale of 
20 fs. If the perfect symmetry of the wave functions is 
even slightly changed due to the excitation of H-H by 
the laser field, the very small spatial separation of the 
electrons leads to their strong Coulombic repulsion all 
over the z-grid. This suggests the reason for the fast 
out-of-phase oscillations of 'atomic' energies E&(t) and 
E-q(£) after the end of the laser pulse [see observation 
(iv) above and Figs. 5(a) and 5(d)]. The fast out-of- 
phase oscillations of the 'atomic' energies occur in the 
case of the entangled initial state and do not occur in 
the case of the unentangled direct-product initial state 
[Figs. 2(a) and 2(e)]. In the case of the direct-product 
initial state, the minimum difference between (zi(t)) and 
(z2(i)) is about 98 a.u. (5.18 nm) [see Figs. 2(c) and 
2(g)], therefore EER is very weak all over the z-grid, ex- 
cept for the special cases of very sharp overlaps of the 
electronic wave functions in the vicinity of protons pa 
and pb (see Figs. 3 and 4). 

A very strong EER taking place in the case of the en- 
tangled initial state should also result in a very efficient 
EER-induced 'sequential' ionization, similar to that de- 
scribed in Sec. Ill for the unentangled direct-product ini- 
tial state of H-H. On the other hand, it was shown above 
that the entanglement of the initial state of H-H does 
not change its ionization probability in comparison to 
the direct-product initial state. Therefore, a closer look 
at the process of EER-induced ionization in the case of 



the entangled initial state is required. To this end, we 
performed several model simulations with the entangled 
initial state excited in the domain of atom A by the laser 
pulses with narrow spatial envelopes defined by Eq. (|16[) 
for various model assumptions of the system-field inter- 
action to be specified below. The results obtained are 
presented in Figs. 7 and 8. 

Electron probabilities P{z\) and P{z2) presented in 
Figs. 7(a) and 7(b) correspond to the end of the 5 fs 
laser pulse with a narrow spatial envelope [Fig. 1(b)] 
which excites both electrons e\ and &i in the vicinity 
of za = —50 a.u. (atom A) and does not affect elec- 
trons ei and ei in the vicinity of zb = 50 a.u. (atom B). 
Apparently, electron probabilities P(z\) and -P(z2), pre- 
sented in Figs. 7(a) and 7(b) at t = 5 fs, are extended 
in comparison to their initial ones in both z\ 2 < Z A 
and 2i 2 > za directions. The laser-induced extension of 
P(zi) and P{z2) into domains zi ; 2 < za gives rise to the 
ionization Ia in the negative directions of the zi^-axes. 
At zi^ < —91 a.u., the wave function is absorbed by the 
imaginary optical potentials. In contrast, at zi,2 > za, 
the laser-driven electrons e\ and e 2 reach the domains 
zi i2 > where they are attracted and accelerated by 
proton pe, localized at zb = 50 a.u., more and more 
efficiently. Due to the electron-proton Coulomb attrac- 
tion, electron probabilities P{z\) and P(z 2 ) increase in 
the vicinity of zip = zb and thus destroy the initially per- 
fect symmetry of the entangled wave function, as clearly 
seen from Figs. 7(a) and 7(b). The extension of elec- 
tron probabilities P{z\) and -P(z2) into the domains of 
#1,2 > zb gives rise to the EER-induced 'sequential' ion- 
ization Ib in the positive directions of the zi^-axes. The 
only problem to be clarified now is EER in the domain 
of atom B, because at a first glance one could conclude 
from Figs. 7(a) and 7(b) that the electrons coming from 
atom A to atom B occupy the domains of z\ t % > zb- 
In order to clarify this issue, we performed two model 
simulations with a 5 fs laser pulse having narrow spatial 
envelope, assuming that only one of the two electrons is 
excited by the laser field in the domain of atom A. The 
results obtained are presented in Figs. 7(c) and 7(d). 

First we assume that only electron e\ is excited by 
the laser field in the domain of atom A. To this end, we 
chose the spatial envelope function F(z) in the interac- 
tion Hamiltonian of Eq. (fj| as follows: F(zi) = Fn(zi) 
and F(z2) = 0, where the narrow spatial envelope F^(z) 
is given by Eq. (fl6|) and illustrated in Fig. 1(b). The 
respective electron probabilities -P(zi) and P{z2) at the 
end of the pulse (t — 5 fs) are shown in Fig. 7(c). Sec- 
ondly, we assume that only electron e 2 is excited by the 
laser field in the domain of atom A. Accordingly, the 
spatial envelope function F{z) in the interaction Hamil- 
tonian of Eq. ([?]) is chosen as follows: F(zi) = 0, and 
F(z 2 ) — F-^(z 2 ). The respective electron probabilities 
P(zi) and P(z 2 ) at the end of the pulse (t — 5 fs) are 
shown in Fig. 7(d). 

It is seen from the results presented in Figs. 7(c) and 
7(d) that the EER-induced ionization of atom B with 
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FIG. 7. Excitation of H-H by the laser pulses with the narrow spatial envelopes from the entangled initial state: electron 
probabilities P (zi) (solid lines) and P(z2) (dashed lines) at the end of the 5 fs laser pulses, (a) and (b) - both electrons ei and 
e2 are excited by the laser pulse in the domain of atom A; (c) - only electron ei is excited by the laser pulse in the domain of 
atom A; (d) - only electron is excited by the laser pulse in the domain of atom A. 



the entangled initial state proceeds similar to that with 
the unentangled direct-product initial state described in 
Sec. Ill above. Specifically, [see Fig. 7(c)] the laser-driven 
electron e\ (coordinate z±) comes from atom A to atom B 
and pushes electron e 2 (coordinate z 2 ) out into the do- 
main z 2 > zb , where the EER- induced ionization of elec- 
tron e 2 takes place in the positive direction of the Z2-axis. 
At the same time, electron e\ which arrived at atom B 
does not affect the probability P(z\) at z\ > 50 a.u.. 
Similarly, [Fig. 7(d)], the laser-driven electron e 2 (coor- 
dinate z 2 ), coming from atom A to atom B, pushes elec- 
tron ei (coordinate z{) into the domain z\ > zb, where 
the EER-induced ionization of electron e\ takes place in 
the positive direction of the zi-axis. Electron ei does not 
affect the probability P(zi) at z^ > 50 a.u.. In a realistic 
case, when both electrons are excited by the laser field 
[Figs. 7(a) and 7(b)], both processes take place simulta- 
neously. 

To identify local maxima of electronic wave functions, 
we next define the probability differences for electrons, 
AP(zi,t) and AP(z2,t), with respect to the initial elec- 
tronic probabilities as follows: 



AP(z u t) = P( Zl ,t) - P(z u t = 0), 
AP(z 2 ,t) = P(z 2 ,t) - P(z 2 ,t = 0). 



(21) 



The probability differences AP(zi) and AP(z2) plotted 
in Fig. 8 are calculated at the end of the 5 fs laser pulse 
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FIG. 8. Excitation of H-H from the entangled initial state 
by the laser pulse with a narrow spatial envelope which af- 
fects only electron ei in the domain of atom A: probability 
differences AP(zi) (solid line) and AP(z2) (dashed line) in 
the domain of atom B at the end of the 5 fs laser pulse. The 
changes of the electron probabilities P(z\) and P(z2) are mea- 
sured with respect to their initial values at t = [see Eq. (|21[) ]. 



with a narrow spatial envelope which excites in the do- 
main of atom A only electron e\ and does not affect 
electron e 2 . It can be concluded from Figs. 7(c) and 
8 that the laser-driven electron e±, coming to the domain 
zi > 0, is attracted and accelerated therein by proton 
Pb and forms at t = 5 fs a local double-peak maximum 
at zb = 50 a.u., where electron e 2 already has its global 
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maximum. Due to a very strong Coulombic repulsion 
of the two electrons in the vicinity of zb = 50 a.u., a 
sharp hole appears in the global maximum of P{z%) at 
the expense of the increasing electron probability P(z2) 
in the domain of z-i > zb, thus giving rise to the EER- 
induced ionization of electron e2 in the positive direction 
of the Z2-axis. A similar process takes place when only 
electron e2 is excited by the laser pulse with a narrow 
spatial envelope, giving rise to the EER-induced ioniza- 
tion of electron e\ in the positive direction of the Zi-axis. 
In a realistic case, when both electrons are excited by the 
laser field in the domain of atom A, both described above 
processes take place simultaneously. 

We conclude from the results presented in Figs. 7 and 
8 that the physical reasons of the energy transfer from 
atom A to atom B and the 'sequential' ionization of 
atom B in the case of the entangled initial state and a 
laser pulse with a narrow spatial envelope, which excites 
only electrons of atom A, are similar to those underlying 
the case of the direct-product initial state: the Coulombic 
attraction of the laser-driven electrons by proton pb , and 
the short-range Coulombic repulsion of the two electrons 
in the vicinity of proton pb, where their wave functions 
strongly overlap. 



V. CONCLUSION 

In the present work we have studied numerically the 
non-Born-Oppenheimer quantum dynamics of two dis- 
tant H atoms (being referred to as atoms A and B) 
with an arbitrary large initial internuclear separation of 
R = 100 a.u. (5.29 nm), which have been excited by 
spatially shaped laser pulses. We have found an effi- 
cient energy transfer from one H atom to the other and 
a 'sequential' ionization of the latter, induced by a short- 
range EER. The short-range EER, taking place e.g. in 
the vicinity of proton pb of atom B, occurs due to the 
preceding long-range LIET from atom A to atom B en- 
hanced by the Coulombic attraction and acceleration of 
the laser-driven electrons by proton pb- 

Both unentangled direct-product atomic states and en- 
tangled molecular states have been used as the initial 
states of the H-H system in our numerical simulations. 
The first case, unentangled atomic states, is more likely 
to arise in an experiment: at a gas pressure of 1 atm., 
for example, the interatomic distance is about 100 a.u.. 
The second case, entangled molecular states, requires an 
additional step: dissociation of the H2 molecule, or gen- 
eration of entanglement between two individual atoms. 
The field strengths of the spatially shaped laser pulses 
were chosen such as not to induce a strong ionization 
and the respective decrease of the overall norm of the 
wave packet of H-H. In stronger fields and with longer 
pulses both the energy transfer and the 'sequential' ion- 
ization can be more efficient than those presented above. 
The parameters of the laser pulses used in the present 
work (in particular the pulse duration, t p = 5 fs, and es- 



pecially the laser carrier frequency, uj = 1.0 a.u.) are not 
yet optimal, and more efforts will be required in order to 
find optimal spatially shaped laser fields suitable for the 
most efficient energy transfer and 'sequential' ionization. 

We have shown that in the case of a narrow spatial 
envelope of the applied laser field, when only electrons 
initially belonging to atom A are excited by the laser 
field, the physical mechanisms of the energy transfer from 
atom A to atom B and the 'sequential' ionization of 
atom B are as follows: (i) the Coulombic attraction of 
the laser-driven electrons by proton pb of atom B, re- 
sulting in the formation of narrow local maxima of the 
electronic wave functions in the vicinity of proton pb, 
and (ii) the short-range Coulomb repulsion of the two 
electrons in the vicinity of proton pb , where their wave 
functions strongly overlap. In a more general case of a 
wide spatial envelope of the laser field, for example a 
broad Gaussian, when both A and B atoms are excited 
by the field simultaneously, the same processes occur in 
the opposite direction as well: the energy is also trans- 
ferred from atom B to atom A and the EER-induced 
'sequential' ionization also occurs in atom A. 

Furthermore, we have shown that long-range entan- 
glement of the initially unentangled direct-product state 
is established by the interaction with a spatially shaped 
laser pulse. Long-range entanglement attracts consid- 
erable interest these days. For example, entanglement 
transfer from dissociated molecules to photons has been 
explored in [l8j], where the concept of transferring the 
quantum state of two dissociated fragments sharing 
internal-translational entanglement to that of two pho- 
tons and vice versa has been put forward. A special kind 
of an entangled state of one electron and two protons in a 
H^~ molecule adiabatically stretched to R = 20 — 100 a.u. 
(1.06-5.29 nm) has been studied recently in [13j. A 
quantum-gate mechanism based on electron spins in cou- 
pled semiconductor quantum dots, separated by 40 nm, 
has been considered in (l9j as an important application of 
long-range entanglement. Such gates provide a source of 
spin entanglement and can be used for quantum comput- 
ers. Experimental generation of entanglement between 
two individual 87 Rb atoms held in two optical tweezers, 
separated by 4 /im, has been reported in (2pj . Finally 
we mention that quantum entanglement among extended 
biomolecules has been explored in [2l| (see also refer- 
ences therein). The results obtained in [21| demonstrate 
that there exists robust quantum entanglement among 
chlorophyll molecules under physiological conditions for 
the case of a single elementary excitation. 

The main characteristic features of the entangled 
molecular initial states in our H-H model, as compared to 
the unentangled atomic direct-product ones, are as fol- 
lows: (i) the immediate response of atom B to the laser 
excitation of atom A, and (ii) the existence of fast out- 
of-phase oscillations of the 'atomic' energies Ep^(t) and 
E#(t) after the end of the applied laser pulse. These 
features should occur in a long-range quantum commu- 
nication among distant quantum systems for the follow- 
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ing reason. The time-dependent 'atomic' energies E^(t) 
and i?B (t) represent the information encoded in the laser 
pulse. Therefore, if only one atom A of the entangled 
A-B pair is excited by the laser pulse, the encoded infor- 
mation reaches atom B (receiver) instantaneously and, 
moreover, can be decoded both during and after the end 
of the pulse acting on atom A. 

We have found that in the case of a narrow spatial 
envelope of an applied laser pulse, which excites only 
electrons belonging to atom A, the 'sequential' EER- 
induced ionization of atom B is more efficient than the 
laser- induced ionization of atom A [see Figs. 2(f) and 
6(e)]. These results suggest important consequences 
for more complex quantum systems, e.g. those com- 
posed of three and more distant sub-systems, such as 
{A„}, n = 1, 2, . . . , N. If we excite by a laser pulse with 
a narrow spatial envelope only sub-system Ai, for ex- 
ample, and if ionization probabilities obey Ia„ < Ia„ + d 
we can stimulate an efficient chain of 'sequential' ioniza- 
tion events with ever increasing probability in the total 
{A„} system. The same may also apply to the energy 
transfer in the {A„} system. To achieve such collective 
energy transfer and 'sequential' ionization, a narrow spa- 



tial shaping of the laser field as reported in Q needs fur- 
ther exploration. To achieve the condition Ia„ < Ia„ + i 
in the {A„} systems with a Gaussian spatial envelope 
of Eq.(Q31), we can always choose zq < such that the 
effective-field amplitudes obey £^ n > £^" +1 . In this case 
we can anticipate that the condition L\„ < Ia„+i is ful- 
filled at the edge of the Gaussian envelope. 

Finally, at a shorter pulse duration and/or at a smaller 
laser carrier frequency, such that the number of optical 
cycles per pulse is N c < 15 (see Ref. @), the role of 
CEP-effects as explored previously in LIET [ll[ offers 
new avenues for energy and information transfer at long- 
range distances. 
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